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A*  general  model  of  charge  transfer  at  the  insulator-semiconduc 
tor  interface  of  metal-insulator-semiconductor  (MIS)  devices  is 
developed.  The  model  is  based  upon  direct  WKB  tunneling  of 
electrons  from  the  semiconductor  substrate  into  localized  trap 
states  within  the  forbidden  gap  of  the  insulator.  Both  an 
arbitrary  distribution  of  the  insulator  states  in  energy  as  well 
as  spatial  variation  in  the  trap  densities  in  the  interface  region 
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are  treated.  Advantage  is  taken  of  the  strong  spatial  decay  of 
the  tunneling  transition  probability  to  obtain  valid  approximate 
solutions  of  the  tunneling  equations  in  a form  directly  amenable 
to  numerical  evaluation.  The  effects  of  internal  electric  fields 
in  the  insulator  are  included  in  a self-consistent  manner.  The 
results  of  several  representative  model  calculations  of 
post- irradiation  annealing  and  field-injection  characteristics 
of  MIS  capacitors  are  described. 
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1 . INTRODUCTION 

A model  is  developed  of  the  charge  transfer  processes  occurring  at 
the  insulator-semiconductor  interface  of  the  gate  insulator  in 
metal-insulator-semiconductor  (MIS)  devices.  This  model  is  based  upon 
direct  WKB  tunneling  of  electrons  into  localized  trap  states  in  the 
insulator.  The  model  is  directed  toward  a general  description  of  both 
annealing  of  holes  that  have  been  trapped  in  an  interface  region  after 
exposure  to  ionizing  radiation  as  well  as  the  bias  instability 
(field -induced  charge  injection)  observed  in  A1203  MIS  devices.  The 
mathematical  framework  of  the  model  is  described  and  representative 
model  calculations  of  postirradiation  annealing  and  field -inject ion 
characteristics  of  MIS  capacitors  are  presented.  A detailed  analysis 
of  rapid  annealing  and  charqe  injection  data  on  Al2 03  MIS  capacitor:- 
in  terms  of  the  tunneling  model  has  been  described  elsewhere.1 


Two  of  the  essential  features*  that  must  be  explained  by  a 
successful  theory  of  charge  transfer  are  (1)  the  strong  bias 
dependence  of  the  field-injection  characteristics  in  the 
AJL2O3  devices,  and  (2)  the  prevailing  logarithmic  dependence  of  the 
transferred  charge  on  time.  The  logarithmic  time  dependence  is  true 
both  in  the  charging  by  field  injection  in  AI2O3  devices  and  in 
annealing  of  positive  charge  after  irradiation  in  both  AI2O3  and 
Sio2  MIS  devices.1'*  The  direct  tunneling  model  does  predict  these 
general  features,  and  it  also  allows  quantitative  information 
concerning  the  oxide  trap  distributions  near  the  interface  'jo  be 
extracted  from  the  analysis  of  experimental  data.  We  should  point  out 
that  the  short-term  annealing  characteristics  of  Si02  MIS  capacitors 

*See  Annotated  Bibliography,  Charge  Injection  and  Post-Irradiation 
Rapid  Annealing  in  NOS  Devices. 

1F.  B.  McLean,  H.  E.  Boesch,  Jr.,  P.  S.  Winokur,  J.  M.  NcGarrity, 
and  R.  B.  Oswald,  Jr.,  IEEE  Trans.  Nucl.  Sci.  NS-21  (1974),  47. 


have  recently  been  shown0-*3  to  be  dominated  by  a stochastic  hopping 
transport  of  holes  to  the  interface,  bat  there  remains  a long-term 
logarithmic  annealing  of  the  po  ^tive  charge  that  occurs  after  most  of 
the  holes  have  transported  to  or  through  the  interface.  This 
long-term  annealing  process  may  be  attributed4  to  relaxation  of  a 
trapped-hole  distribution  in  the  interface  region  via  tunneling  of 
elections  from  the  Si  substrate  with  subsequent  recombination  with  the 
trapped  holes. 

The  modeling  of  the  charge  transfer  at  the  gate  insulator 
interfaces  employs  single  electron  WKB  theory 5 for  electrons 
tunneling  directly  from  the  electronic  energy  bands  of  the  substrate 
material  (normally  Si)  or  gate  metal  into  localized  electron  states 
with  energies  in  the  forbidden  gap  of  the  insulator.  It  draws  to  some 
extent  upon  the  idea  behind  a certain  class  of 
metal -nitride-oxide-semiconductor  (MNOS)  memory  transistors.*  The 
physical  basis  of  the  model  is  also  similar  to  that  underlying  the 
theory  of  the  frequency  response  of  surface  states  in  MIS  systems. ^ 
The  theory  is  developed  formally  in  a general  way  and  includes  in 
principle  both  energy  conserving  and  energy  nonconserving  processes 
(e.g.,  phonon-assisted  tunneling)  as  well  as  three-dimensional 
effects.  The  localized  insulator  states  are  either  electron  traps,  as 
for  field  injection,  or  occupied  hole  traps  in  postirradiation 
annealing.  The  model  treats  both  an  arbitrary  energy  distribution  of 

*See  Annotated  Bibliography,  Theory  of  MNOS  Memory  Transistors . 

t See  Annotated  Bibliography,  Theory  of  Frequency  Reponse  of  Surface 
States. 

2F.  B.  McLean,  G.  4.  Ausvean,  Jr.,  H.  E.  Boesch,  Jr.,  and 
J.  M.  McGarrity,  J.  Appl.  Phys.  49  (1976),  1529. 

3Jf.  E.  Boesch,  Jr.,  F.  B.  McLean,  J.  M.  McGarrity,  and 
G.  A.  Ausman,  Jr.,  IEEE  Trans.  Nucl.  Scl.  NS- 2 2 (1975),  2163. 

J.  Meier,  Summaries  of  Papers,  IEEE  Annual  Conf.  on  Nucl.  and 
Space  Radiation  Effects,  Seattle  (1972),  215. 

5C.  B.  Duke,  Tunneling  in  Solids,  Academic  Press,  New  York  (1969), 
p.  30. 
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insulator  states  as  well  as  spatial  variation  in  the  trap  densities  in 
the  interface  regions.  The  known  energy-band  structure  of  silicon  (or 
the  gate  metals)  convoluted  with  the  capture  cross  sections  of  the 
insulator  states  is  used  as  the  "supply  function"  for  the  tunneling 
electrons.  The  effects  of  internal  electric  fields  in  the  insulator 
are  included  in  a self-consistent  fashion  by  solving  Poisson's 
equation  along  with  the  tunneling  equations  (i.e . , the  change  in  the 
internal  electric  fields  due  to  the  accumulation  of  tunneling  charge 
is  accounted  for) . 

In  section  2 the  general  theory  is  formulated  and  illustrated  by 
the  exact  solution  for  a simple  case,  namely,  that  of  a single  trap 
level  with  uniform  spatial  distribution.  In  section  3 the  general 
theory  is  developed  analyticall ■'  to  the  point  where  the  solutions  can 
be  obtained  numerically  by  simple  quadrature.  This  is  accomplished  by 
means  of  a set  of  simplifying  approximations  which  are  justified  on 
physical  grounds.  In  section  4 details  of  the  numerical  procedure  are 
described,  and  model  calculations  for  several  representative  cases 
important  in  practice  are  presented.  First,  the  field-injection 
characteristics  for  AI2O3  devices  are  described  for  two 
cases:  (1)  where  a thin  Si02  layer  is  present  between  the  Si 
substrate  and  the  deposited  AI2O3,  and  (2)  where  no  Si02  layer  is 
present.  The  important  qualitative  features  distinguishing  these  two 
cases  are  pointed  out.  Second,  the  charge  relaxation  characteristics 
following  irradiation  due  to  tunneling/recombination  of  a trapped- hole 
distribution  in  the  interface  region  is  discussed.  Particular 
emphasis  is  placed  upon  the  effects  produced  by  spatial  variation  of 
the  trapped-hole  distribution.  In  the  last  section  the  main  results 
are  summarised,  and  a brief  discussion  is  given  of  the  application  of 
the  model  to  the  analysis  of  experimental  data. 
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2.  FORMAL  THEORY 


Let  N(E,x)  be  the  arbitrary  distribution  of  insulator  traps,  both 
in  energy  E and  distance  x from  the  tunneling  interface 
(insulator-semiconductor  interface) . Energy  values  will  be  referenced 
to  the  valence  band  edge  of  the  insulator  (fig.  1).  Let  n(E,x,t)  be 
the  distribution  of  traps  occupied  at  time  t.  We  assume  that  at  time 
t = 0 all  traps  are  unoccupied,  i.e.,  n(E,x,0)  a 0.  Next,  we  define  a 
function  q(E,x,t)  that  gives  the  probability  per  unit  time  that  an 
electron  will  make  a tunneling  transition  from  the  substrate  material 
to  a single  unoccupied  trap  of  enerqy  E located  at  the  distance  x at 
time  t.  The  function  g(E,x,t)  is  essentially  a Green's  function  for 
the  tunneling  transition  rate;  it  is  time  dependent  when  the  internal 
electric  field  in  the  interface  region  is  time  dependent,  as  is  the 
case  when  the  charge  distributions  are  changing  significantly  due  to 
the  charge  transfer  process.  The  tunneling  transition  rate  T(E,x,t) 
to  traps  of  energy  E at  distance  x from  the  interface  is  simi  ly  the 
product  of  q(E,x,t)  and  the  unoccupied  trap  density 


T(E,x,t)  - g(E,x,t) lN(E,x)  - n(E,x,t) ] 


But  this  is  also  equal  to  the  time  rate  of  change  of  the  occupied  trap 


distribution 


dn(E,x,t) 


T(E,x,t) 


the  solution  of  which  is 


n(E,x,t)  ■ N(£,x)  Jl  - exp 


[C at' 


g(E,x,t' 


The  total  tunneling  current  per  unit  area  J(t)  is  then  given  by  the 
double  integral  of  T(E,x,t)  over  all  E and  x.  Using  equation  (3),  we 
have  the  formal  solution  for  the  tunneling  current 
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INSULATOR  SEMICONDUCTOR 

Figure  1.  Schematic  band  diagram  of  the  insulator-semiconductor 
interface  in  tin  MIS  system  at  zero  bias.  Distance  is 
measured  with  respect  to  the  interface  and  energies  are 
measured  relative  to  E^,  the  valence  band  edge  of  the 
insulator.  Et  represents  a localized  trap  state  in  the 
insulator,  which  acts  as  a recipient  for  electrons 
tunneling  from  the  semiconductor  substrate. 


J(t) 


E 

sj  dx/  c dEg(E,x,t)N(E,x)3xppy^t  dt’g(E,x,t' )J 


(4) 


where  e is  the  electronic  charge  and  E*  is  the  conduction  band  edge  of 
the  insulator  (measured  relative  to  the  insulator  valence  band  edge). 
The  total  charge  transferred  at  time  t is  simply 

Q(t)  - ft  dt'J(t')  (5) 
o 


It  should  be  noted  that  J(t)  is  not  the  current  that  would  be  measured 
by  an  external  circuit,  but  is  the  charge  flux  at  the  tunnelinq 
interface.  The  externally  measured  current  would  be  much  smaller,  due 
to  the  small  displacement  of  the  tunneling  charge  relative  to  the 
thickness  of  the  insulator  film. 
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The  tunneling  equations  must  he  solved  self -consistently  with 
Poisson's  equation  to  account  for  the  relaxation  of  the  energy  bands 
due  to  the  tunneled  charge.  This  band  relaxation  is  described  by  the 
time  dependence  of  the  electrostatic  potential  $ (x,t)  given  by 

4>(x,t)  = 4k(x)  - ejx  £(x',t)dx'  (6) 

o 

where  is  the  image  potential  (important  only  very  near  the 

interface)  and  £(x,t)  is  the  x component  of  the  electric 
field.  £(x,t)  is  found  by  the  solution  to  Poisson's  equation 

£(x,t)  - £ (t)  + [*  dx'  P (7) 

° Jo  H 

whe re  £ it)  is  the  field  at  x * 0,  e.  is  the  dielectric  constant  of 

O 1 

the  insulator,  and  p(x,t)  is  the  charge  distribution  in  the  insulator 
given  by 

E1 

p(x,t)  * P0(x)  + e|  c dE  n(E,x,t)  (8) 

where  p (x)  is  the  initial  charge  distribution. 


Equations  (4)  and  (5)  solved  self-consistently  with  Poisson's 
equation  thus  constitute  the  formal  solution  to  the  tunneling  problem. 
In  essence  the  problem  is  reduced  to  determining  the  Green's  function 
g(E,x,t).  Various  approaches  are  possible  to  determine  g(E,x,t)  with 
varying  degrees  of  sophistication,  such  as  whether  one  considers  one- 
or  three-dimensional  tunneling,  elastic  (energy  conserving)  or 
inelastic  tunneling.  However,  the  only  viable  starting  point  for  any 
approach  is  the  one-electron  NKB  approximation.  The  essential  results 
of  this  approximation  are  that  we  consider  only  one-electron 
transitions  from  the  Bioch  band  states  of  the  substrate  solid  to  the 
localised  insulator  states  and  that  the  probability  of  making  a 
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transition  to  a state  at  distance  x from  the  interface  is  proportional 
tc  the  exponential  factor  exp(-Z (E,x,t) ) where  Z is  essentially  a 
scaled  distance  given  by 


Z(E,x,t)  = 2 J X dx'  B(E,x't) 
o 


B(E,x,t) 


(9) 

(10) 


Here,  V_(x,t)  is  the  x-  and  t-dependent  potential  barrier  through 
B * 

which  the  electrons  tunnel  and  m is  the  effective  mass  of  the 

* * 

electrons  in  the  barrier  region.  (In  general,  m^  can  also  be  energy 
dependent.)  Thus,  within  the  WKB  approximation  g(E,x,t)  has  the 
general  form 

g(E,x,t)  ■ aiE,x,t)e  z^E'x*t^  (11) 


where  everything  but  the  exponential  decay  factor  is  lumped  into  a 
"prefactor"  a(E,x,t).  The  exponential  term  may  loosely  be  thought  of 
as  controlling  the  penetration  of  the  electrons  into  the  insulator 
and  a as  being  the  "supply  function"  for  the  tunneling  electrons, 
i.e.,  the  number  of  electrons  per  unit  time  attempting  the  transition. 
The  reason  for  writing  g in  the  form  of  equation  (11)  is  that  except 
for  band-edge  effects — to  be  discussed  in  more  detail  later— the 
prefactor  term  a(E,x,t)  is  a relatively  slowly  varying  function  of 
both  E and  x as  compared  with  the  exponential.  As  we  see  below,  this 
fact  can  be  used  to  great  advantage  in  effecting  the  integration  of 
equation  (4) . 

1 

{ 

\ 

In  order  to  obtain  a qualitative  feelinq  for  the  nature  of  the 
solution,  consider  a very  simple  case  for  which  J(t)  and  Q(t)  can  be 
directly  evaluated  without  making  further  approximations.  Namely,  we 
consider  the  case  of  a single  trap  level  of  energy  E with  uniform 


11 


spatial  distribution  N.  Only  one-dimensional  tunneling  is  considered 
and  the  effects  of  the  internal  electric  fields  are  neglected  (valid 
only  for  small  charge  densities  at  low  applied  bias  levels). 
Then  a is  a constant  and  Z = 2f5x  is  a function  only  of  x with 


Direct  integration  of  equations  (4)  and  (5)  yields 

J,tl  ' (l  ‘ e‘at) 

Q(t)  - [ei  (at)  + in  (at)  + y] 


(12) 

(13) 


where  y = 0.57721  ...  is  Euler's  constant  and  Ej  is  the  exponential 
integral  defined  fcy6 

00  -f  • 

El  (T)  - / ~~  df  (14) 

T 

The  limiting  forms  of  Ej(t)  are 


for  t < 1 


for  t > 1 


(15) 


Note  in  particular  that  Ej(t)  goes  to  zero  rapidly  for  t > l. 


6 A#.  A bramowitz  and  I.  A.  Stegun,  Handbook  of  Mathematical 

Functions,  Dover  Publications,  Inc.,  New  York  (1965). 
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Some  of  the  main  features  of  the  results,  equations  (12)  and 
(13),  are  as  follows:  The  time  scale  of  the  tunneling  current  is  set 

by  the  value  of  a.  At  early  times,  t <<  a-1#  J(t)  and  Q(t)  have  the 
simple  form 


J(t) 


(16a) 


for  t « 1/a 


Q(t) 


- Nat  /.  at  \ 

26  \ 4/ 


(16b) 


Thus,  J (t)  -*■  constant  and  Q(t)  increases  linearly  with  t as  t-*- 0.  At 
t ~ a-1  the  form  of  J (t)  changes  from  that  of  equation  (16)  to  an 
inverse  time  relationship,  and  for  > a-1  we  have  sinply 


J(t) 


20t 


Q(t)  s Jjj-  ttnat  + y) 


(17a) 

for  t >>  1/a 

(17b) 


In  this  time  regime  the  characteristic  logarithmic  dependence  of  the 
total  transferred  charge  is  obtained.  If  3 is  known  then  a plot  of  Q 
versus  log  time  yields  directly  the  trap  density. 


The  solution  to  the  simple  case  under  consideration  can  be 
extended  immediately  to  an  arbitrary  trap  distribution  in  energy  N{E), 
but  where  the  distribution  still  is  spatially  uniform.  Still  only 
one-dimensional  tunneling  is  considered,  and  band  bending  due  to 
internal  electric  fields  is  neglected.  Here  a and  6 become  functions 


« 


of  E,  but  there  is  no  coupling  between  different  values  of  E,  so  that 
equations  (12)  and  (13)  become  simple  integrals  over  the  energy 
distribution 


yE1  N(E)  f,  _-a(E)tj 

J<t)  ■/  c “ Ll  ' e J 


O 

,i 


Q(t)  = / C dE  [ei  (a(E)t)  + in(a(E)t)  + y] 


28  (E) 


(18) 


(19) 


Note  that  for  t » l/a(E)  for  all  E values  of  interest,  the 
characteristic  t”1  dependence  of  J(t)  and  ln(t)  dependence  of  Q(t)  is 
still  obtained. 


Before  considering  the  integration  of  equations  (4)  and  (5)  for 
J(t)  and  Q(t)  in  the  general  case,  we  examine  the  etructure  of  the 
occupied  trap  distribution  n(x,t)  and  the  transition  rate  function 
T(x, t)  for  the  simple  one-level  case.  This  offers  valuable  insight  on 


how  to  proceed  with  the  general  case.  Substituting  g ■ ae 
equations  (3)  and  (1)  yields 


-8x 


n(x,t) 


T(x,t) 


l - .-'‘I 


-28x  -ate 
Nae  e 


-28x 


into 


(20) 


(21) 


In  figure  2 the  normalized  distribution  n/N  versus  x is  plotted  at 

values  of  time  spaced  at  order-of -magnitude  intervals.  Time  is  scaled 

in  units  of  a-1  and  x is  scaled  in  units  of  (28) The  dashed  curve 

simply  reflects  the  exponential  decay  of  the  tunneling  probability 

with  x (not  to  scale) . The  important  point  of  figure  2 is  that  at  a 

given  time  the  trap  distribution  is  essentially  saturated  out  to  the 

vicinity  of  some  point  x (to  be  rigorously  defined  below)  about  which 

m 

the  distribution  falls  nearly  to  zero  in  a distance  of  the  order  of 
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Figure  2.  The  normalized  occupied  trap  distribution  as  a function  of 
distance  from  the  tunneling  interface  at  decade  increments 
in  time  for  a single  trap  level.  Time  is  scaled  in  units 
of  a-1  and  x in  units  of  (2B)_1.  The  dashed  curve  reflects 
the  exponential  decay  of  the  relative  tunneling  probability 
with  distance. 

several  scaled  units.  The  point  x for  t > a-*  obviously  increases 

m 

with  time  in  a logarithmic  fashion.  Since  the  transition  rate 

function  T(x,t)  is  proportioned  to  the  product  of  the  number  of 

unoccupied  traps  and  the  exponential  decay  factor,  we  see  that  T(x,t) 

will  be  peaked  rather  sharply  in  the  vicinitv  of  the  off  of  the 

occupied  trap  distribution.  In  fact,  we  shall  define  x (t)  to  be  the 

hi 

point  at  which  the  tunneling  transition  rate  is  maximized  at  time  t; 
i.e,,  we  find  x for  which  3T(x,t)/3x  » 0,  with  the  result 


xa(t) 


0 for  t < a"’ 

^ in  (at)  for  t > a-1 


(22) 
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This  again  describes  the  logarithmic  time  dependence  of  the  total 
transferred  charge.  In  fact,  comparing  equations  (17)  and  (22)  for 
t >>  a-1  we  see  that  Q(t)  can  be  written 


«<«  ■ “[«„<'>  + 2«] 


(23) 


The  term  Y/2  3,  which  is  of  the  order  of  0.5  a for  typical  values  cf  3, 

accounts  for  the  asymmetry  of  n(x,t)  around  x (t) . For  t » a 1 the 

m 

distribution  n(x,t)  has  a universal  shape  relative  to  the  point  x (t) 

m 

and  the  same  is  true  for  the  profile  (but  not  magnitude)  of  T(x,t). 

Figure  3 shows  the  normalized  distribution  plotted  versus  x-x  (t) , and 

m 

also  shows  the  profile  of  the  transition  rate  about  x (t) . We  see 

m 

that  T (x,t)  is  indeed  peaked  rather  sharply  about  x (t)  with  a 

m 

width  at  half -maximum  (FWHM)  of  three  scaled  units  (23)_1 . For  the 
situations  to  be  considered  in  the  typical  MIS  structures,  3 generally 
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DISTANCE  X-X.  MOV  POINT  OF  MAXIMUM  TUNNELING  RATE  (UMTS  Of  (l A) 

Figure  3.  The  normalized  occupied  trap  distribution  and  relative  tun- 
neling rate  plotted  against  distance  from  the  point  cf  the 
maximum  tunneling  rate. 
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lies  in  the  range  0.5  to  1.0  inverse  angstrom  units.  Hence,  the  full 

width  at  half  maximum  (FWHM)  of  T(x,t)  is  of  the  order  1.5  to  3.0  A. 

Thus,  at  any  time  t,  most  of  the  tunneling  "action"  is  taking  place 

to  localized  states  lying  in  a narrow  region  of  width  several 

angstroms  centered  around  the  distance  x (t)  from  the  interface. 

m 

3.  APPROXIMATE  SOLUTION  IN  TOE  GENERAL  CASE 


The  formalism  of  the  previous  section  will  now  be  applied  to  the 

solution  of  the  tunneling  problem  in  the  general  case, 

♦ 

namely,  where  we  want  to  consider  (1)  the  effects  of  internal 
electric  fields,  (2)  possible  spatial  variations  in  the  trap 
densities,  (3)  three-dimensional  tunneling,  find  (4)  realistic  forms 
for  the  prefactor  term  a(E,x,t)  of  equation  (11)  as  calculated  within 
some  model.  In  the  general  case  it  is  not  possible  to  obtain  am 
analytic  solution.  Equations  (4)  and  (5)  must  be  solved  numerically 
along  with  Poisson's  equation,  and  a(E,x,t)  in  equation  (11)  is 
calculated  within  some  reasonable  approximation.  However,  the 
logarithmic  dependence  of  the  transferred  charge  on  time  is 
characteristic  of  direct  tunneling  into  localized  states  and  carries 
over  generally,  but  with  some  modifications,  to  more  complicated 
cases.  Modifications  to  the  ln(t)  behavior  are  due  mainly 
to:  (1)  the  relaxation  of  the  energy  bands  in  the  insulator  due  to 

the  tunneled  charge,  which  is  particularly  important  for  tunneling 
from  states  near  a band  edge  in  the  semiconductor,  and  (2)  spatial 
variations  in  the  trap  density  in  the  interface  region.  In  this  last 
case  it  is  possible  to  map  out  the  actual  spatial  distribution  of  the 
trap  density  near  the  interface  by  analysis  of  the  slopes  of  annealing 
or  injection  curves  (as  long  as  the  band  relaxation  effects  are 


The  essential  point  is  that  even  in  the  general  case,  as  long  as 

the  electrostatic  potential,  ij)(x,t),  and  N(E,x)  vary  sufficiently 

slowly,  it  is  the  exponential  term  in  equation  (11)  that  dominates  the 

spatial  dependence  of  g(E,x,t)  and,  consequently,  of  T(E,x,t). 

Therefore,  the  qualitative  picture  depicted  in  figure  3 is  still 

valid.  In  particular,  at  time  t the  tunneling  transition  rate  for 

electrons  of  energy  E is  peaked  at  a distance  x (E,t)  with  a spread 

n 


around  x (E,t)  of  the  order  of  several  angstroms, 
in 


We  may  take 


advantage  of  this  feature  by  expanding  all  terms  in  the  integral  of 

“Z 

equation  (4)  except  the  exponential  terms,  e , about  the  point 

x (E,t)  (keeping  in  mind  the  energy  dependence  of  x ),  which  is  found 
m m 

by  maximizing  T(E,x«>t)  with  respect  to  x.  As  long  as  <Hx,t)  and 

O 

N(E,x,t)  vary  slowly  over  distances  of  the  order  of  ~3  A,  it  is 

sufficient  to  retain  only  the  first  or  second  terms  of  the  expansions, 

and  the  evaluations  of  J(t)  and  Q(t)  are  simplified  considerably. 

* 

Omitting  the  details  of  the  calculations  here,  we  give  the  essential 
results  for  the  case  where  only  the  leading  order  terms  in  the 
expansions  are  retained 
following  conditions. 


The  point  x (E,t)  is  found  to  satisfy  the 
m 


If  x (E)  denotes  the  point  that  maximizes 
m 

T(E,x,t)  at  t * 0,  then  for  t < t (E),  x (E,t)  ■ x°(E),  where  the  time 

o m m 

t (E)  is  determined  by  the  solution  of  the  equation 
o 


l 

/ 


dt’g(E,x°(E),t') 


(24) 


For  t > t (E) , x (E, t)  is  found  to  satisfy  the  identity 

o m 


r 


dt'g(E,xn(E,t)  ,t^  - 1 


(25) 
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Differentiating  equation  (25)  with  respect  to  t yields  a differential 

equation  for  x (E,t) 
in 


3x  (E,t) 
m 


3t  26q(E) 


hs)  9(E,xm(E.t),t) 


where 


B <E) 
o 


•Sn-fl' 


Equation  (26)  can  be  simply  integrated  forward  in  time  to  give  x (E,t) 

m 


i rt 

x (E,t)  = x°(E)  + / dt'g(E,x  (E,f),f) 

m m 2B  (E)  \ m / 


o'  ' t (E) 
o 


It  is  important  to  keep  track  of  the  time  arguments  in  equations  (24) 

to  (28) . We  point  out  that  x°(E)  may  be  nonzero  in  either  of  two 

in 

cases:  (1)  In  a double  insulator  system  where  there  are  no  traps  in 

the  first  insulator — that  nearest  to  the  semiconductor  substrate  (see 
below  for  further  discussion) , and  (2)  for  energy  states  near  a band 
edge  of  the  substrate  solid,  where  the  states  at  the  surface  may  be  in 
the  band  gap  of  the  semiconductor,  but  due  to  band  bending  by 
electrostatic  fields,  tunneling  may  be  allowed  to  states  deeper  in  the 
insulator . 


Once  x (E,t)  is  determined,  the  tunneling  current  J(t) 

m 

and  the  total  transferred  charge  Q(t)  are  found  by  carrying  out  Taylor 

expansions  about  x (E,t)  and  performing  the  integration  over  x in 

m 

equation  (4)  and  over  t'  in  equation  (5)  analytically.  The  results  in 
the  general  case  can  be  expressed  in  a form  analogous  to 
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equations  (12)  and  (13)  for  the  simple  case  considered  earlier  (one 
energy,  one -dimensional  tunneling,  no  band  bending).  We  define  a 
scaled  time  parameter  i(E,t). 


x (E , t)  ' h(E,t)exp 


2&0(E)  f “ 

J om 

m 


x (E,t) 


dxN(E,x)! 


(29) 


where 


h(E,t)  =/*  dt'g(E,xm(E,t)  ,t') 

A ' 


(30) 


Note  h(E,t)  = 1 for  t > t (E)  by  the  identity  (eq  (25)).  In  terms  of 

o 

the  scaled  time  x(E,t),  J(t)  and  Q(t)  can  be  expressed  as 

J(t)  . . L ^'V‘)  h . ,-TB.tO  (31) 

Jtt)  ej  aE  2go(E)  h(E, t)  L1  e (31) 


Q(t) 


/N(E,x ) 

* 2g(E)/  CEl(T)  + tn(T)  + *] 

o 


where  x » x (E,t)  and  t » t (E , t ) are  understood, 
m m 


(32) 


Of  course,  the  tunneling  equations  (24)  to  (32)  are  to  be  solved 
simultaneously  with  Poisson's  equation.  however,  we  again  take 
advantage  of  the  nature  of  the  occupied  trap  distribution  as  shown  in 
figure  2.  Namely,  it  is  a good  approximation  to  assume  that  at  time  t 
traps  of  energy  E are  essentially  saturated  to  the  point  xn(E,t)  and 
empty  beyond.  Then  an  analytical  solution  to  4(x,t)  in  terms  of  a 
simple  quadrature  over  E can  be  effected.  Let  4q(x)  be  the 
electrostatic  potential  at  t ■ 0,  which  is  determined  by  the  applied 
bias,  work  function  difference,  image  potential,  and  the  initial 


?0 


charge  distribution  in  the  insulator.  We  assume  that  the  change  in 
the  potential  with  time  is  due  solely  to  the  tunneled  charge.  Then, 
at  time  t, 

4>(x,t)  * 4>0(x)  “ \ Q(t)  J7  * hJ dEN(E»t^F(X}Xn(B)  »xm(B#t))  (33* 


where 


x (E, t) 


N(E,t) 


Jr*  iis 
m 

X°(E) 

n 


dxN(E,x) 


(34) 


and 


for  x < Xi 
for  xi  < x < X2 


f(xix1(x2)-{(1/2)(x  - xi)2 

(*2  “ xi)  (x l— — for  x > x2 


(35) 


with  x2>xj.  Also  in  equation  (33),  ■ e^/l^  i»  the  insulator 
geometrical  capacitance,  where  is  the  insulator  dielectric 
constant,  and  is  the  insulator  thickness. 

By  taking  advantage  of  the  peaked  structure  of  T(E,x,t)  around 

x (E,t) , we  have  been  able  to  transform  the  general  solution  to  the 
in 

tunneling  problem,  equations  (4)  and  (5) , into  a form  directly 
amenable  to  numerical  calculation.  One  simply  integrates 
equation  (28)  forward  in  time  in  conjunction  with  simple  quadratures 
over  energy.  All  that  remains  is  to  determine  an  appropriate  form  for 
the  Green's  function,  g(E,x,t),  or  more  to  the  point,  determine  an 
appropriate  form  for  the  prefactor  a(E,x,t)  in  equation  (11).  Ne  have 
not  carried  out  a complete  quantum-mechanical  calculation,  but  we  have 
calculated  g within  a semiclassical  treatment  that  contains  all  the 
salient  features  of  a more  rigorous  treatment.  In  particular,  our 
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treatment  includes  the  density  of  states  of  the  semiconductor 

(therefore  band -edge  effects  are  included)  and  accounts  for 

three-dimensional  effects.  We  take  the  electrons  within  a band  in  the 

semiconductor  to  be  simulated  by  a free  electron  gas  with  an  effective 

mass  m*  appropriate  to  the  band  electrons.  We  treat  the  electrons 

«-  £ 

striking  the  boundary  as  particles  having  a WKB  probability  factor  e 
for  penetration  into  the  barrier  region  x > 0,  which  is  weighted 
according  to  the  normal  component  of  the  electron  velocity  v^  in  the 
manner 

<(E,x,t)  * ~ Z‘,£,x,t)  (36) 

x 

where  Z(E,x,t)  is  given  by  equation  (9)  and  v is  the  total  electron 
velocity.  This  weighting  procedure  simulates  the  more  rigorous,  but 
much  more  complicated,  quantum  mechanical  result5  in  which  electrons 
striking  the  surface  at  near-normal  incidence  have  a much  larger 
penetration  probability  to  a distance  x than  those  with  an  appreciable 
velocity  component  parallel  to  the  surface.  For  electrons  of  energy  E 
the  calculation  of  g involves  a summation  over  all  possible  wave 
vectors  having  a positive  x component.  We  give  the  final  result  for 
elastic  tunneling  of  electrons  from  the  semiconductor  valence  band, 
taken  to  be  isotropic,  to  a localized  trap  of  energy  E located  at  a 
distance  x from  the  interface 

* 

g(E,x,t)  - a(E)0(E,)E’F(Z)e"Z  (37) 


5C.  fl.  Duke,  Tunneling  in  Solids,  Academic  Press,  New  York  (1969), 
p.  30. 
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where 


E*  * Ey  - E + *(x,t) 


0(E' ) 


[l 


for  E'  > 0 
for  E'  <0 


Z “ 2x 


* ) 

p 


F(Z)  = (1/2)  [l  - Z + eZE!(Z)] 


(38) 


(39) 


(40) 


(41) 


and  o(E)  is  an  effective  capture  cross  section  per  unit  energy  for  the 
traps  of  energy  E and  has  to  be  determined  by  comparison  with 

experimental  data.  The  m*  appearing  in  equation  (37)  is  appropriate 

* 

for  the  semiconductor  valence  band  and  the  m in  equation  (40)  is  that 
appropriate  for  electrons  in  the  potential  barrier  region.  For 
multiple  overlapping  valence  bands--as  in  the  case  for  Si— the  m*  in 
equation  (37)  is  simply  replaced  by  a sum  over  the  m*  for  each  band. 
The  function  F(Z)  accounts  for  the  three-dimensional  effects;  it  has 
the  limiting  forms 


F(Z)  * 


1/2 

for  Z -►  0 

(1/Z)  (] 

-§) 

for  Z » 1 

(42) 


For  trap  states  near  the  interface,  essentially  all  electrons  in  the 
hemisphere  for  v^ > 0 contribute  equally  to  the  tunneling  probability, 
whereas  for  states  far  removed  from  the  interface,  only  those 
electrons  striking  the  surface  with  near-normal  incidence  contribute 
significantly  to  the  tunneling.  Thx ' decrease  in  the  number  of 
effective  electrons  is  reflected  in  tl , Z 1 dependence  of  F(Z)  for 
large  Z.  In  writing  equation  (40)  we  have  assumed  that  4(x,t)  can  be 
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approximated  by  a straight  line  in  the  insulator.  However,  this 
linear  approximation — being  made  within  the  square  root — introduces 
little  error  as  long  as  $(x,t)  is  reasonably  smooth  and  the  total 
change  in  $ in  the  tunneling  region  is  not  comparable  to  the 
magnitude  of  E*  - E,  which  is  assumed  to  be  true  for  our  present 
applications. 

The  result  for  g(E,x,t)  has  been  given  for  electrons  tunneling 
from  the  valence  band  of  the  semiconductor.  Analogous  results  hold 
for  electrons  tunneling  from  the  conduction  band  as  well  as  for 
tunneling  at  the  metal  insulator  surface.  The  results  can  also  be 

extended  to  include  tunneling  from  an  anisotropic  band  described  by  an 

* 

effective  mass  m^  for  electron  wave  vectors  normal  to  the  interface 

* 

and  by  an  effective  mass  mu  for  electron  motion  parallel  to  the 

interface.  The  essential  result  for  the  case  Z >>  1 is  that  the  ra*  of 

* * * , 

equation  (37)  is  replaced  by  ®i.(®ii/®].)  • For  the  conduction  band  of 
silicon  oriented  in  the  100  direction  this  leads  to  a reduction  in  the 
tunneling  rate  by  a factor  of  •*100. 

The  theory  has  been  developed  explicitly  for  situations  for  only  a 
single  insulator.  But  the  theory  is  readily  extended  to  the  double 
insulator  system,  as  is  thought  to  be  the  situation  for  AI2O3  HIS 
devices  where  a thin  intermediary  SiOx  layer  of  width  10  to  30  A lies 
between  the  Si  substrate  and  the  AI2O3  insulator  film.  We  assume  that 
trapping  states  in  the  intermediary  film  are  unimportant  and  that 
tunneling  occurs  from  the  band  states  of  the  semiconductor  through  the 
intermediary  film  and  into  localized  states  in  the  second  insulator 
film.  If  x is  the  width  of  the  intermediary  film  and  tfte  coordinate 
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system  is  located  so  that  x = 0 refers  to  the  interface  between  the 
two  insulators,  so  that  N(E,x)  ^0  only  for  x > 0,  then  equation  (37) 
becomes 


* 

g(E,x,t)  - o(E)0(E')E'F(z  + Z \e_  (2  + Zo) 

(2»Ti) 3 ' °> 


where 


- E + <Mx,t) 


_ »(o,t)^ 


H 


(43) 


(44) 


In  2 , m*  and  E1  refer  to  the  tunneling  effective  mass,  and 
o t c 

conduction  band  edge,  respectively,  of  the  intermediary  insulator, 
whereas  in  equation  (40)  for  2 they  refer  to  those  quantities  in  the 
second  insulator.  Also,  in  expression  (33)  for  the  electrostatic 
potential,  a term 


r eolt) 

e 1 o 


x 

o 


must  be  added  to  account  for  the  potential  drop  across  the 

intermediary  film.  Here,  £ (t)  is  the  electric  field  just  inside  the 

o 

second  insulator  and  the  ratio  e^/c j of  the  dielectric  constants  of 
the  two  insula*  urs  accounts  for  the  discontinuity  at  the  boundary 
x = 0 between  the  two  (different)  dielectric  media. 


4.  MODEL  CALCULATIONS 

In  this  section  model  calculations  will  be  presented  based  on  the 
numerical  solution  of  the  tunneling  equations  developed  in  the 
previous  section.  These  calculations  ar  ! for  relatively  simple  model 
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systems  in  order  to  exhibit  the  major  features  and  predictions  of  the 
theory.  A detailed  comparison  of  the  model  with  charge  injection  data 
on  A1203  MIS  capacitors  has  been  described  elsewhere.1 


The  model  systems  used  in  the  calculations  were  chosen  to  be 
representative  of  several  cases  important  in  practice.  First,  the 
field  injection  characteristics  of  Al203  MIS  capacitors  are  calculated 
for  two  cases:  one  in  which  a thin  Si02  layer  is  present  between  the 
Si  substrate  and  the  deposited  A1203  film,  and  the  other  in  which 
there  is  no  Si02  layer  present.  The  energy  band  diagrams  for  these 
two  cases  are  depicted  schematically  in  figure  4 for  zero  and 

positive  applied  bias  conditions.  The  energy  positioning  of  the 
0 b e 
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Figure  4.  Schematic  energy  band  diagrams  for  model  calculations  (top 
diagrams:  flat  band  condition;  bottom  diagrams:  positive 

applied  bias) ; (a)  single-insulator  system,  (b)  double- 
insulator system — localized  traps  associated  with  second 
(extended)  insulator,  and  (c)  double-insulator  system — 
traps  associated  with  intermediate  insulator  or  with  inter- 
face region  between  two  insulators. 


1F.  B.  McLean,  H.  E.  Boesch,  Jr.,  P.  S.  Winokur,  J.  M.  McGarrlty, 
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electron ' trap  band  (to  be  discussed  in  more  detail  below),  is  such 
that  at  flat-band  conditions  the  electron  traps  have  energies  lying 
opposite  to  the  forbidden  gap  of  the  silicon  substrate  (top  diagrams 
in  fig.  4) . When  positive  bias  is  applied  to  the  MIS  structure,  these 
states  shift  down  in  energy  relative  to  the  Si  bands,  and  for 
sufficiently  large  bias,  direct  tunneling  will  occur  into  them  from 
the  Si  valence  band  (bottom  diagrams  of  fig.  4).  Finally,  the  charge 
relaxation  characteristics  following  irradiation  are  calculated  due  to 
tunneling  of  electrons  into  the  insulator  and  recombining  with  a 
trapped  hole  distribution  in  the  interface  region.  This  case  is 
also  depicted  schematically  in  figure  4.  The  hole  traps  are  assumed  to 
have  energies  which  at  ilat-band  conditions  are  below  the  Si  valence 
band  edge.  Though  normally  unoccupied,  these  states  become  populated 
by  capturing  holes  generated  by  ionizing  radiation.  Holes  trapped 
sufficiently  close  to  the  interface  subsequently  can  serve  as  receptor 
states  for  electrons  tunneling  from  the  Si  valence  band.  This 
tunneling/recombination  process  may  be  an  important  mechanism  for  the 
annealing  of  radiation-induced  positive  charge  at  the  Si  interface  in 
both  AI2O3  and  Si02  MIS  devices.  In  the  calculations  presented  here, 
emphasis  is  placed  upon  the  effects  produced  by  spatial  variation  of 
the  trapped  hole  distribution.  In  certain  instances  it  may  be 
possible  to  profile  the  spatial  distribution  of  the  trap  density  near 
the  interface  by  appropriate  analysis  of  experimental  data. 

The  primary  input  to  the  numerical  calculations  is  the  spatial  and 
energy  distribution  of  the  insulator  trap  states.  For  the 
calculations  reported  here,  the  energy  and  spatial  distributions  are 
assumed  to  be  separable,  i.e.,  N(E,X)  * N(E)f(x).  The  energy 
distribution  N (E)  is  taken  to  be  gaussian,  and  the  spatial 
distribution  f(x)  to  be  an  analytic  form  that  is  integrable  in 
expressions  such  as  equations  (29)  and  (34).  The  calculations  proceed 
by  the  integration  of  equation  (28)  for  X (E,t)  forward  in  time  along 
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with  the  integral  of  Poisson's  equation,  equation  (33),  and  utilizing 

equations  (37)  to  (41)  describing  the  form  of  the  Green's  function 

g(E,x,t)  employed  here.  Once  X (E,t)  is  known,  equation  (29)  is 

m 

evaluated  for  the  scaled  time  parameter  t(E,t),  along  with  the 
function  h(E,t).  Then  the  tunneling  current  J(t)  and  accumulated 
charge  Q(t)  are  calculated  by  simple  integrations  over  the  trap  energy 
distribution  with  equations  (31)  and  (32) . For  the  double  insulator 
case  the  modifications  described  at  the  end  of  section  4 cure  used,  we 
now  discuss  the  results  of  the  representative  model  calculations. 

4.1  Field-Injection  Characteristics 

We  first  present  the  model  calculations  of  field- inject ion 
characteristics  for  several  cases  appropriate  to  AI2O3  MIS  systems, 
assuming  that  direct  tunneling  is  the  mechanism  at  work.  As  mentioned 
earlier,  for  the  modeling  we  consider  a band  of  electron  traps  which 
at  flat-band  conditions  lies  adjacent  to  the  forbidden  gap  of  the  Si 
substrate.  This  configuration  was  chosen  on  the  basis  of  photo- 
depopulation data  presented  by  Harari  and  Royce7  which  indicated  that 
such  a band  (their  so-called  "c"  band)  is  a major  recipient  of 
field- injected  electrons.  Further,  photodepopulation  studies  carried 
out  on  samples  used  in  a previous  HDL  study1  indicate  that  at  room 
temperature  and  for  bias  levels  up  to  about  25  v above  the  flat-band 
voltage,  almost  all  the  field-injected  electrons  ended  up  in  states  in 
this  same  energy  region.  In  general,  other  electron  trap  levels  may 
come  into  play  and  perhaps  other  processes,  particularly  at  higher 
bias  levels  where  tunneling  from  the  Si  conduction  band  into  higher 
trap  levels,  or  into  the  insulator  conduction  band  is  possible. 
Nevertheless,  the  essential  features  of  the  field-injection  charac— 

1F.  B.  McLean , H.  E.  Boesch,  Jr.,  P.  S.  Winokur,  J.  M.  McGarrity, 
and  R.  B.  Oswald,  Jr.,  IEEE  Trans.  Nucl.  Sci.  NS-21  (1974),  47. 

1E.  Harari  and  B t S.  H.  Royce , IEEE  Trans.  Nucl.  Sci.  NS-20 
(1973),  280. 
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teristics  can  be  extracted  based  upon  the  following  simple  idea.  This 
band  of  electron  traps  is  nearly  empty  at  flat-band  conditions  due  to 
the  unavailability  of  electrons  for  direct  tunneling  into  them,  but  as 
positive  bias  is  applied  the  band  bends  downward  with  respect  to  the 
Si  band  structure  (see  fig.  4a  and  4b).  At  sufficiently  high  bias, 
direct  tunneling  is  possible  from  the  Si  valence  band  into  the 
localized  electron  states  in  the  insulator.  This  idea  forms  the  basis 
for  the  theory  behind  a certain  class  of  MNOS  memory  transistors. * 

In  particular,  for  our  model  calculations  we  consider  a trap 
band  of  uniform  spatial  density  of  5.0  x 1019  cm-3 (based  on 
experimental  data1) , but  having  a gaussian  energy  distribution  peaked 
at  flat -band  condition  at  an  enerqy  of  0.6  eV  above  the  Si  valence 
band  edge  and  with  a FWHM  of  0.2  eV.7  The  Si  valence  band  edge  is 
located  4.1  eV  below  the  AI2O3  insulator  conduction  band.  The  total 

O 

oxide  thickness  is  taken  to  be  1000  A,  and  the  effective  mass  of  the 
tunneling  electrons  is  taken  as  0.4  times  the  free  electron  mass.8 

Figure  5 presents  the  fie Id- injection  curves  appropriate  for 

P 

a double  insulator  situation  where  a thin  Si02  film  20  A thick  is 
postulated  to  exist  between  the  Si  substrate  and  the  extended 
AI2O3  insulator.  This  is  thought  to  be  the  case  by  many 


*See  Annotated  Bibliography , Theory  of  MNOS  Memory  Transistors. 
lF.  B.  McLean,  H.  E.  Boesch,  Jr.,  P.  S.  Winokur,  J.  M.  McGarrity, 
and  R.  B.  Oswald,  Jr.,  IEEE  Trans „ Nucl.  Sci.  NS-21  (1974),  47. 

7E.  Harari  and  B.  S.  H.  Roy ce,  IEEE  Trans.  Nucl.  Sci.  NS- 20 
(1973),  280. 

SM.  Lenslinger  and  E.  H.  Snow,  J.  Appl . Phys.  40  (1967),  278. 
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TIME  (S) 

Figure  5.  Model  calculations  of  field  injection,  AVFB  versus  time,  for 
double-insulator  MIS  system  with  20-A  Si02  film  between  Si 
substrate  and  extended  AI2O3  insulator.  Electron  trap  dens- 
ity ■ 5.0  * 10 19  cm”3,  electron  trap  band  width  - 0.2  eV, 
and  mean  energy  is  0.6  eV  above  Si  valence  band  (applied 
bias  for  each  curve  shown  in  parentheses) . 


investigators7  when  AI2O3  is  directly  deposited  onto  a Si  substrate. 
Here  we  assume  that  the  electron  traps  are  associated  with  the 
A1203  insulator  material.  Figure  5 shows  the  shift  in  flat-band 
voltage  (AV __(t)«Q(t))  versus  time  after  application  of  positive  bias 

FB 

for  a series  of  bias  voltages  (indicated  by  the  numbers  in 
parentheses).  Noticeable  injection  begins  at  a bias  level  of  about 
8 V relative  to  flat  band,  and  the  injection  rapidly  picks  up  with 


7ff.  tfarari  and  8.  5.  H.  Roy ce,  IEEE  Trans.  Nucl . Sci.  NS-20 
(1973),  280. 
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increasing  bias.  We  note  that  the  injection  is  effectively  controlled 
by  the  intermediate  Si02  film;  it  is  the  potential  drop  across  this 
film  which  sets  the  injection  threshold.  Further,  the  time  scale  for 
tunneling  is  controlled  by  its  thickness;  as  a consequence  the 
injection  curves  for  all  bias  levels  begin  (on  a log  scale)  at  about 
the  same  time.  The  other  general  feature  to  note  here  is  that  the 
injection  curves  initially  break  into  the  characteristic  ln(t) 
behavior,  but  as  time  progresses  they  tend  to  flatten  out  (especially 
evident  at  the  higher  bias  levels) . This  apparent  saturation  effect 
is  a direct  consequence  of  the  relaxation  of  the  bands  due  to  the 
tunneled  charge,  and  this  relaxation  is  amplified  across  the  thin 
Si(>2  film  due  to  the  difference  in  dielectric  constants  of  A1203  and 
Si02«  (We  take  the  relative  dielectric  constant  to  be  3.85  for 
SiC>2  and  8.5  for  AI2O3.) 

Figure  6 shows  the  injection  curves  for  a single  insulator 
situation;  a similar  trap  structure  as  in  the  previous  case  is 
considered,  but  the  thin  SiC>2  film  is  omitted.  Now,  within  our  model 
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Figure  6.  Calculated  field-injection  characteristics  for  single- 
insulator MIS  system.  The  same  electron  trap  parameters 
were  used  as  those  in  figure  5. 


the  initial  tunneling  rate  is  determined  by  the  position  of  the  closest 
available  states  for  direct  tunneling  from  the  Si  valence  band,  and  for 
this  case  this  distance  is  inversely  proportional  to  the  electric  field 
at  the  interface.  As  the  field  is  increased  the  insulator  energy  bands 
bend  downward  more  sharply  and  this  distance  decreases  (see  fig.  4a), 
leading  to  an  increased  initial  tunneling  rate.  Consequently,  as  the 
bias  is  increased  the  injection  curves  begin  breaking  into  the  ln(t) 
behavior  at  increasingly  earlier  times.  Other  points  to  note  are  that 
the  slopes  of  the  curves  in  the  log  regime  are  nearly  parallel  (but  with 
some  fanning  out  evident)  and  that  the  relaxation  effect  due  to  the 
tunneled  charge  is  considerably  reduced  because  of  the  absence  of  the 
intermediate  Si02  film  of  lower  dielectric  constant. 

Figure  7 shows  the  calculated  bias-stress  curves  for  the  single 
and  double  insulator  cases  where  the  flat-band  shift  10^  s after 
application  of  the  bias  is  plotted  versus  the  applied  bias.  Both 
curves  exhibit  a pronounced  threshold  for  injection,  as  is  observed 
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Figure  7.  The  calculated  flat  band  voltage  shift  after  application  of 
bias  for  1000  s plotted  versus  applied  gate  bias  for  the 
single-  and  double- insulator  model  systems. 
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experimentally.  The  threshold  voltage  bias  is  considerably  smaller  for 
the  double  insulator  system  (for  the  model  parameters  used  in  the 
present  calculations)  than  for  the  single  } insulator.  This  can  be 
directly  attributed  to  the  enhanced  potential  drop  (or  band  bending) 
across  the  intermediate  film  as  a result  of  the  mismatch  in  dielectric 
constants,  thereby  producing  smaller  values  of  X° (E)  in  the  double 
insulator  system.  This  would  no  longer  be  the  case  at  sufficiently  high 

O o 

bias  levels  where  X < 20  A for  the  single  insulator  system,  but  this 

m 

occurs  outside  the  bias  range  considered  here.  The  other  major 
qualitative  feature  distinguishing  the  two  cases  is  that  there  is  a 
"slowing  down"  of  the  field  injection  at  the  higher  bias  levels  for  the 
double  insulator  case  as  compared  with  the  bias-stress  curve  for  the 
single  insulator,  where  strong  injection  continues  with  increasing  bias. 
The  difference  again  can  be  attributed  to  the  more  pronounced  relaxation 
effect  in  the  double  insulator  system. 

Next,  we  discuss  the  effect  of  the  trap  bandwidth  on  the 
injection  characteristics  for  the  single  insulator.  In  figure  6,  where 
the  bandwidth  is  taken  to  be  0.2  eV,  the  transition  region  from  the 
initial  flat  slope  (on  a log  plot)  to  the  logarithmic  regime  is 
relatively  narrow.  However,  if  the  bandwdith  is  increased  to  0.6  eV,  as 
for  the  injection  characteristics  shown  in  figure  8,  this  transition 
region  is  broadened  considerably.  This  broadening  is  due  basically  to 
the  increased  spread  in  the  closest  states  available  for  direct 
tunneling  between  the  top  and  bottom  of  the  electron  trap  band.  In 
practice,  by  appropriate  analysis  of  the  injection  curves  for  this  case, 
information  can  be  obtained  concerning  the  energy  structure  of  the  trap 
band.  To  see  this,  refer  to  figure  9,  where  at  flat-band  conditions  the 
bottom  and  top  of  the  trap  band  are  taken  to  be  Aj  and  A2,  respectively, 
relative  to  the  Si  valence  band  edge.  The  transition  time  ti  is  defined 
as  the  time  at  which  a noticeable  flat-band  shift  due  to  the  injection 
process  is'  first  obtained;  it  corresponds  approximately  to  the  time  at 
which  a significant  quantity  of  electrons  have  tunneled  into  the  nearest 


Figure  8.  Calculated  field-injection  characteristics  for  a single- 
insulator MIS  system  illustrating  the  effect  of  increased 
electron  trap  band  width  (0.6  eV) . 

ANALYSIS  OF  INJECTION  CURVE:  SINGLE  INSULATOR 


Figure  9.  Determination  of  electron  trap  band  parameters  from  a typical 
experimental  injection-dominated  transient  flat  band  voltage 
response  curve. 


available  states  for  tunneling  (at  xj).  These  states  are  associated 
with  the  lower  edge  of  the  trap  band  (at  A}).  The  second  transition 
time,  t2»  is  defined  as  the  time  at  which  the  injection  curve  approaches 
true  logarithmic  behavior.  It  corresponds  to  the  time  at  which  the 
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closest  available  states  (at  X2)  or  traps  near  the  upper  edge  of  the 

trap  band  (at  A2)  have  become  saturated,  and  after  which  the  tunneling 

transition  rate  to  all  energy  levels  in  the  band  decays  as  t-1  (and 

hence  AV  increases  as  In  (t».  Now,  within  the  WKB  approximation  the 
FB 

characteristic  tunneling  time  t.  to  states  located  A.  above  E at 

1 IV 

flat-band  conditions  is  given  by  t^  = t?  exp(26xj  where  B is  an 

averaged  value  of  the  WKB  integrand  (eq  (9) ) over  the  energy  range  of 

concern,  and  x^  is  the  distance  required  for  the  electrostatic  potential 

to  change  by  an  amount  A^  from  its  value  at  the  interface.  Since 

x.  = A ./£  where  £ is  the  electric  field  in  the  interface  region,  the 
110  o 

expression  for  t^  can  be  rewritten  as 

t. 

In  ^5  = 2BA./<?o  . (45) 

Hence,  an  estimate  of  the  values  of  Aj  and  A2  can  be  obtained  by  plotting 
the  log  of  tj  and  t2  versus  £~ 1 . The  intercepts  of  such  plots  at 

£q  1 * 0 can  be  used  to  obtain  estimates  of  the  capture  cross  sections  of 

the  traps.  1 

To  conclude  this  section  we  note  the  injection  characteristics 
computed  by  using  the  direct  tunneling  model  are  qualitatively  different 
for  the  single  and  double  insulator  systems,  and  in  practice  it  should 
be  readily  apparent  which  case  one  is  observing.  However,  if  one  is 
dealing  with  a double  insulator  system  vhere  the  electron  trap?.,  instead 
of  being  associated  with  the  extended  insulator  as  in  figure  5,  are  in 
fact  associated  with  the  intermediate  film  (see  fig.  4c) , then  the 
qualitative  features  of  the  injection  characteristics  would  correspond 
to  that  of  the  single  insulator  situation,  as  shown  in  figures  6 and  8. 
A more  subtle  analysis  is  then  required  to  distinguish  between  the 
cases. 

1F.  B.  McLean,  H.  E.  Boesch,  Jr.,  P.  s.  winokur,  J.  M.  McGarrity,  and 
R.  B.  Oswald,  Jr.,  IEEE  Trans.  Nucl.  Sci.  NS-21  (1974),  47. 


4.2  Post-Irradiation  Charge  Relaxation 


Here  are  discussed  the  model  calculations  of  rapid  annealing 
via  tunneling/recombination  into  a single  band  of  trapped  holes  near  the 
Si  interface  following  a pulse  of  ionizing  radiation.  The  energy  levels 
of  the  trapped  hole  states  are  taken  to  be  adjacent  to  the  Si  valence 
band.  In  particular,  the  peak  of  the  band  is  taken  to  lie  1.8  eV  below 
the  Si  valence  band  edge  with  a FWHM  of  0.6  eV,  which  are  reasonable 
values  based  on  photodepopulation  studies.7 *'9  We  restrict  our  attention 

O 

to  the  single  insulator  case  (of  thickness  1000  A) , since  the  annealing 
characteristics  are  qualitatively  similar  vftiere  a thin  intermediate 
insulator  film  is  present,  the  only  important  difference  being  a shift 
in  the  onset  time  for  significant  tunneling.  We  use  values  for  the 

tunneling  effective  mass  and  dielectric  constant  appropriate  to  that  of 

* 

Si02  (mt  » 0.4  me,  = 3.85).  It  is  assumed  that  after  the  MIS 
capacitor  has  been  subjected  to  the  pulse  of  ionizing  radiation  the  hole 
trap  band  is  saturated  in  the  vicinity  of  the  interface  to  a distance 
beyond  tunneling  range  with  a resulting  flat-band  shift  of  -10.0  V.  The 
focus  of  attention  here  is  the  manner  in  which  spatial  variations  in  the 
trapped  hole  distribution  affect  the  annealing  characteristics.  To 
illustrate  the  effect  we  assume  a simple  exponential  form  for  the 
spatial  dependence  of  the  distribution  near  the  interfaces 

N (x)  - NQeXX  (46) 

where  is  the  density  of  traps  at  the  interface  (integrated  over  all 
energies).  We  take  Nq  - 2 * 1019  cm”3,  a typical  value  estimated  from 
experimental  data. 


7E.  Harari  and  B . S.  H.  Royce,  IEEE  Trans.  Nucl.  Sci.  US-20  (1973), 
280. 

Harari,  S.  Wang  and  B.  S.  H.  Royce,  J.  Appl.  Phys.,  (1975), 

1310. 
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Figure  10  shows  the  annealing  curves  for  several  values 

of  X from  0.1  to  -0.2  in  units  of  20  (0  being  evaluated  at  the  peak  of 

the  trap  band).  Clearly,  the  annealing  reflects  directly  the  spatial 

variation  of  the  hole  distribution.  This  is  a direct  consequence  of  the 

nature  of  the  tunneling/recombination  process , namely,  that  there  is  a 

distinct  tunneling  "front"  x (t),  which  moves  into  the  insulator  as  log 

m 

time,  with  the  states  behind  the  front  filled  and  those  beyond  it  still 

empty.  At  a time  t the  tunneling  takes  place  predominantly  to  traps 

located  within  a few  angstroms  of  x (t),  and  hence  reflects  the  density 

m 

of  traps  in  that  region.  (Recall  fig.  2 and  3.)  More  specifically,  it 
is  the  slope  of  the  annealing  curve  vfoich  is  directly  proportional  to 
the  density  of  traps.  Hence,  for  a constant  spatial  density  (X  « 0.0) 
the  annealing  curve  is  very  nearly  linear  with  log (t)  for  t > a"1.  For 
nonzero  values  of  X the  initial  annealing  behavior  is  very  similar  to 


Figure  10.  Calculated  annealing  curves,  showing  effect  of  spatial 
variation  in  traps  distributions.  is  in  units  of 

a-1  and  distance  in  units  of  (20) “1. 


the  constant  density  case,  but  with  increasing  time  the  slopes  of  the 
curves  increase  or  decrease  reflecting  either  an  increasing  or 
decreasing  trap  density  away  from  the  interface. 

Therefore,  if  it  is  certain  that  the  tunneling/recombination 

mecnanism  is  responsible  for  charge  annealing,  the  spatial  profile  of 

the  trapped  hole  distribution  near  the  interface  can  be  obtained  by 

simple  analysis  of  the  experimental  annealing  curves,  provided  that  the 

variation  of  the  tunneling  rate  with  energy  over  the  trap  bandwidth  is 

sufficiently  small.  This  is  the  situation  for  the  present  model 

calculations  for  which  the  energy  levels  of  the  hole  traps  are  located 

well  into  the  Si  valence  band,  and  therefore  the  exponential  WKB  spatial 

decay  factor  dominates  the  transition  rate.  If  the  trap  band  were 

located  near  a band  edge  of  the  substrate,  this  would  no  longer  be  true 

as  there  is  then  a strong  variation  with  energy  in  the  density  of 

electrons  available  for  tunneling.  The  tunneling  "front"  would  then  be 

smeared  out  over  a considerable  range  due  to  the  energy  dependence  of 

X (t) , and  the  annealing  curves  would  no  longer  be  simply  related  to  the 
m 

spatial  density  of  traps.  By  the  same  token  there  is  little  bias 
dependence  in  the  annealing  characteristics  shown  in  figure  10  for  the 
present  model  system.  However,  if  the  trap  band  were  near  a Si  band 
edge,  a strong  bias  dependence  of  the  annealing  would  exist  similar  to 
the  strong  bias  dependence  of  the  field- injection  characteristics 
discussed  earlier. 

5.  SUMMARY 

A general  model  of  charge  transfer  at  the  insulator-semiconductor 
interface  has  been  described.  The  model  is  based  upon  direct  WKB 
tunneling  of  electrons  from  the  semiconductor  substrate  into  localized 
trap  states  in  the  insulator.  Advantage  was  taken  of  the  strong  spatial 
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decay  of  the  tunneling  transition  probability  to  obtain  valid 
approximate  solutions  of  the  tunneling  equations.  These  solutions  are 
of  a form  directly  amenable  to  numerical  evaluation  in  cases  where  there 
is  a distribution  of  the  localized  insulator  states  both  in  energy  and 
in  space.  The  effects  of  internal  electric  fields  in  the  insulator  cure 

included  self-consistently  by  solving  Poisson's  equation  along  with  the 
tunneling  equations.  The  results  of  several  model  calculations 
representative  of  systems  important  in  practice  were  described.  Some  of 
the  important  points  discussed  are: 

(a)  The  direct  tunneling  of  electrons  into  trap  states  in  the 
insulator  yields  the  generally  observed  log  time  feature  of  charge 
transfer  characteristics. 

(b)  The  slopes  of  the  charge  transfer  curves  in  the  logarithmic 
regime  are  proportional  to  the  spatial  density  of  the  traps. 

(c)  Spatial  variations  in  the  trap  densities  modify  the  log(t) 
dependence. 

(d)  The  field-injection  characteristics  cure  qualitatively  different 
for  single  and  double  insulator  systems. 

(e)  A detailed  analysis  of  the  field-injection  curves  can  yield 
information  on  the  energy  distribution  of  the  traps. 

In  a recent  publication,  the  direct  tunneling  model  ms  applied  in 
some  detail  to  the  analysis  of  field-injection  data  in  AI2O3  MIS 
capacitors. 1 It  was  shown  that  electron  injection  under  positive  bias 


1F.  B.  McLean , H.  E.  Boesch,  Jr.,  P.  S.  winokur,  J.  M.  McGarrity,  and 
R.  B.  Oswald,  Jr.,  IEEE  Trans.  Nucl.  Sci.  NS-21  (1974),  47. 
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(as  well  as  trapping  of  radiation-generated  holes)  is  dominated  by  an 
interface  transition  region  at  the  Si-Al203  interface.  The  extent  of 
this  region  from  the  Si  substrate  as  inferred  from  the  data  is  of  the 
order  of  25  K.  The  field  injection  characteristics  were  qualitatively 
explained  by  direct  tunneling  from  the  Si  valence  band  into  electron 
traps  in  a single  insulator.  Quantitative  consistency  between  the 
single  insulator  model  and  the  results  of  photodepopulation  measurements 
after  field  injection  was  obtained  if  it  was  assumed  that  the  dielectric 
properties  of  the  transition  region  were  similar  to  that  of  the  SiC>2 , 
and  further  that  the  electron  traps  were  associated  with  the  transition 
region  rather  than  being  situated  in  the  AI2O3.  Thus,  the  application 
of  the  model  in  this  case  led  to  an  improved  understanding  of  the  nature 
of  the  interface  region  and  of  the  factors  invovled  in  the  field- 
injection  process. 

A detailed  application  of  the  direct  tunneling  model  to  charge 
annealing  after  irradiation  has  not  yet  been  undertaken.  There  is 
recent  experimental  evidence2'3  in  Si02  MIS  capacitors  that  indicates 
that  the  dominant  factor  in  the  charge  relaxation  of  that  system  is  a 
slow,  temperature-activated  transport  of  the  radiation-generated  holes 
through  the  insulator  to  the  Si  interface  where  most  of  the  holes 
(-90  per-ent)  pass  on  to  the  substrate.  The  remaining  fraction  of 
holes,  however,  appears  to  be  trapped  near  the  interface.  A 
tunneling/recorabination  mechanism  involving  these  trapped  holes  may  be 
responsible  for  the  long-term  annealing  process  observed  to  occur  in 
these  samples  for  hours  or  even  days  after  the  irradiation.  It  is 
planned  to  investigate  the  long-term  annealing  in  some  depth. 


2F.  B.  McLean,  G.  A.  Ausman,  Jr.,  H.  E.  Boesch,  Jr.,  and 
J.  M.  McGarrity,  J.  Appl.  Phys.  £9  (1976),  1529. 

E.  Boesch,  Jr.,  F.  B.  McLean,  J.  M.  McGarrity,  and 
G.  A . Ausman,  Jr.,  IEEE  Trans.  Nucl.  Sci.  NS- 2 2 (1975),  2163. 
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